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We present theoretical results on the high-temperature phase stability and phonon spectra of paramagnetic 
bcc iron which explicitly take into account many-body effects. Several peculiarities, including a pronounced 
softening of the [110] transverse (T,) mode and a dynamical instability of the bcc lattice in harmonic 
approximation are identified. We relate these features to the a-to-y and y-to-6 phase transformations in 
iron. The high-temperature bcc phase is found to be highly anharmonic and appears to be stabiUzed by the 
lattice entropy. 

The extraordinary magnetic and metallurgical properties of elemental iron (Fe) have attracted the attention of 
mankind for several thousand years already. In particular, iron is the main ingredient of steel which is one of 
the most important construction materials of modern civilization'. In this respect we still live in the "iron 
age". Iron is also of primary interest for geophysics as the main constituent of the Earth's core^"^. Furthermore, 
iron is perhaps the most famous example of allotropy in metals. Indeed, iron is a very complicated material where 
magnetic order""" and correlation effects'^"^" play an important role. For these reasons the electronic and lattice 
properties of iron are still far from understood. 

Iron has at least four allotropes: the a, y, and 5 phases at low, and the 8 phase at high pressure, respectively. It is 
particularly remarkable that the low-temperature structure of iron, the ferromagnetic a. phase, is not close-packed 
but body-centered cubic (bcc). This already indicates that iron has highly unusual features. In fact, among the 
elements iron is characterized by the unique property of having two stable bcc phases: the aforementioned a 
phase, which is stable below 1185 K, and the 5 phase which is stable above 1670 K up to the melting curve. 
Sandwiched in between those two phases is the y phase ("austenite") which has a closed-packed face-centered 
cubic (fee) structure. It is understood that magnetism (iron is ferromagnetic below the Curie temperature T^ ~ 
1043 K) stabilizes the low- temperature a phase by compensating for the too low cohesive energy. Indeed, the low- 
temperature bcc phase does not exist for nonmagnetic Ru and Os, which lie below Fe in the periodic table. 
Moreover, iron makes a transition to a hexagonal closed-packed (hep) structure (e phase) once ferromagnetism is 
suppressed by pressure. Clearly it is the interplay between magnetism of the electrons and the lattice degrees of 
freedom which determines the complexity of the low-temperature phase diagram of iron. 

At high temperatures, in the absence of ferromagnetic long-range order, the situation is even less clear. In general, 
polymorphic metals which have both bcc and closed-packed structures (e.g., Fe, Ca, Sr, Ti, Zr, Ce, Pu) are known to 
melt fi-om the bcc phase. As argued by Zener''', this can be caused by a somewhat lower Debye temperature in a bcc 
lattice. Namely, the bcc phase has a low-energy Tj phonon [transverse (110)] branch, whose frequencies vanish in the 
case of pair-wise, nearest-neighbor interatomic interactions. Such a soft-mode behavior is associated with a large 
amplitude of vibrations. It leads to a high lattice entropy, which lowers the free energy and thereby stabilizes the bcc 
lattice. This behavior is clearly revealed in prototypical systems with metaUic chemical bonds, such as alkali and 
alkaline-earth metals^'*, as well as in the high temperature bcc phase of the group 3 (Sc, Y, La) and 4 (Ti, Zr, Hf) 
metals^^. By contrast, the phonon modes in bcc transition metals like Cr, Mo, and W, where the valence d-band is 
almost half-filled and the chemical bond is covalent, are hard. Moreover, these metals are not polymorphic. For a 
general discussion of the effects of electronic structure on lattice properties of metals and alloys see Ref 26. 
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It is quite surprising to note that in the case of iron the experi- 
mental information about electronic and dynamical lattice prop- 
erties at high temperatures is still very limited. In fact, to the best 
of our knowledge results on the lattice dynamical properties of 5 iron 
are presently not available. In particular, soft-mode behavior of the 
phonon spectra has neither been reported at the a-to-8 phase trans- 
ition, nor in the y phase^^ ™. Thus the lattice dynamics, in particular 
its evolution at high temperatures and near the a-to-y and y-to-S 
phase transitions in iron, remain unexplored. 

An overall explanation of the electronic and lattice properties of 
iron at high temperatures requires a formalism which can treat both 
paramagnetic and magnetically ordered states on equal footing and 
can describe local moments above Tg. Therefore conventional band- 
structure techniques cannot be used to describe the a-to-y phase 
transition in paramagnetic iron, although they provide a good quant- 
itative description of the properties of ferromagnetic a iron' By 
contrast, the structural transition actually occurs at a temperature 
Tstruc which lies about 150 K above T^. Moreover, these techniques 
find the bcc phase to be elasticaUy unstable in the absence of mag- 
netization. Today we know^"'^' that an improved description of these 
transitions can only by obtained when the existence of local moments 
above Tc is taken into account. 

To clarify these fundamental questions we performed microscopic 
calculations of the structural phase stability and the phonon spectra 
of iron within the so-called LDA + DMFT approach''" It merges 
ab initio techniques for the calculation of band structures, such as the 
local density approximation (LDA) or the generalized gradient 
approximation (GGA), with dynamical mean-field theory 
(DMFT)"" - a state-of-the-art computational scheme for the cal- 
culation of the electronic structure of strongly correlated materials. 
This combined approach allows one to compute the properties of 
correlated materials in both their paramagnetic or magnetically 
ordered states. Several applications of this scheme have provided 
very good quantitative descriptions of the electronic, magnetic, and 
structural properties of iron" In particular, the LDA + DMFT 
method implemented with plane-wave pseudopotentials^''' has 
recently been employed to compute the equilibrium crystal structure 
and phase stability of iron near the a-to-y phase transition^"'^'. 
Electronic correlations were found to be essential for the explanation 
of the electronic and structural properties of paramagnetic iron. For 
example, the a-to-y phase transition was found to occur at —1.3 Tc, 
i.e., well above the magnetic transition, in agreement with experi- 
ment. Moreover, the calculated phonon spectra of iron near the a-to- 
y phase transition were also shown to be in good agreement with 
experiment. 

Results 

We computed the temperature evolution of the lattice dynamical 
properties of iron near the a-to-y and y-to-5 phase transformations 
within LDA + DMFT for Coulomb interaction parameters U = 
1.& eV and / = 0.9 eV as obtained by earlier theoretical and experi- 
mental estimations'^'^*. Our calculations of the Curie temperature in 
bcc Fe yield Tc = 1650 K, which is in reasonable agreement with the 
experimental value 1043 K given the local nature of the DMFT 
approach"*^ and the approximate density-density form of the 
Coulomb interaction^^. To extract the phonon frequencies, the 
dynamical matrix was evaluated using the method of frozen pho- 
nons" implemented within LDA + DMFT. By introducing small 
displacements of atoms from their equilibrium positions, we deter- 
mined the phonon frequencies at given q-points. The phonon dis- 
persion curves were obtained by using a standard Born- von Karman 
(BvK) force constant model with interactions expanded up to the 
fifth-nearest neighbor. In the following all results will be presented as 
a function of the reduced temperature T/Tc, where Tc is the calcu- 
lated Curie temperature (1650 K). 
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Figure 1 | Calculated phonon dispersion curves for bcc iron near the 
Curie temperature Tc. The results are compared with neutron inelastic 
scattering measurements at 1043 K. 

Our results for the phonon dispersion curves of iron near the 
magnetic phase transition temperature Tc are shown in Fig. 1. The 
calculated phonon dispersions exhibit the typical behavior of a bcc 
metal. The phonon frequencies are overall positive, both in the fer- 
romagnetic and paramagnetic phase. This implies the stability of the 
bcc lattice structure of iron in accordance with experiment. We also 
note the remarkably good quantitative agreement between the cal- 
culated phonon dispersions of bcc iron and the experimental data^" 
taken at T ~ Tc (a detailed comparison at 1.2 Tc was reported 
earlier^'). In the given temperature range the calculated phonon dis- 
persions show only a rather weak temperature-dependence. 
However, there is an anomaly in the transverse (110) acoustic mode 
(Tj mode) along the T-N branch near the bcc-fcc phase transition 
temperature at T ~ 1.2 Tc. This can be ascribed to a dynamical 
precursor effect of the bcc-to-fcc phase transition which occurs at 
Tstruct ~ 1-3 Tc in our calculations™. 

Next we discuss the lattice dynamics of paramagnetic bcc iron in 
more detail. Our results (Fig. 2), obtained in the harmonic approxi- 
mation, clearly indicate that at T ~ 1.4 Tc, i.e., above the a-to-y 
phase transition temperature, the bcc phase is dynamically unstable. 
The origin of the instability lies in the Tj mode in the T-N direction 
whose frequency becomes imaginary near the N-point. At the same 
time other phonon branches remain stable and are weakly temper- 
ature dependent near the transition point. At even higher tempera- 
tures, T ~ 1.8 Tc, in the 8 phase, the bcc lattice structure is 
dynamically unstable in the whole F-N branch, again due to the Tj 
mode, with an additional anomaly near the (2/3,2/3,2/3) point. The 
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Figure 2 | Calculated phonon dispersions of paramagnetic bcc iron near 
the a-to-y and y-to-6 phase transitions for different temperatures. 
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Figure 3 | Calculated total energy of paramagnetic iron as a function of displacement (in au) for the longitudinal (L) and the two transverse ( Tj and Tj) 
phonon modes at the N point for the bcc lattice, and the transverse (111) mode at the L-point for the fee lattice, respectively. 



latter anomaly may be interpreted as a precursor effect similar to the 
one occurring at the P-co transition in Ti and Zr^^"*". Therefore our 
calculations reveal the crucial importance of the Tj phonon mode for 
the explanation of the structural phase stability of paramagnetic bcc 
iron, in agreement with the Zener's hypothesis^'. Anharmonic lattice 
effects must be included to obtain a complete description of the 
lattice dynamics of 5 iron. This situation is quite similar to that in 
the 8 phase of Pu, which has a bcc structure and which is dynamically 
unstable in the harmonic approximation'". It also seems to be the 
case in other bcc transition metals, such as Sc, Y, La, Ti, Zr, Hf, and 
'pjj25,42 HQ^gyg]-^ jj^ tjjg g phase of Fe and the s phase of Pu the 
situation is very special since, both, anharmonic lattice and elec- 
tron-electron correlation effects are highly relevant, making the 
problem particularly difficult to investigate. 

To illustrate the importance of anharmonic effects in bcc iron, we 
calculated the total energy profiles for the longitudinal (L) and the 
two transverse (Tj and T2) modes at the N-point. The results are 
shown in Fig. 3, as a function of the corresponding atomic displace- 
ment for three different temperatures. We compare our results with 
those obtained for the lowest phonon, the transverse (111) mode, of 
the fee phase. It is seen that the two bcc phonon modes L and T2 are 
rigid and do not change across the ot-to-y phase transition. Similarly, 
the fee (111) phonon mode is stable and is almost temperature inde- 
pendent. By contrast, the Tj mode is strongly affected by temper- 
ature: while it is stable at T ~ 1 .2 Tc < Tstruo i-e-i below the calculated 
a-to-y phase transition temperature, a shallow double-well potential 
develops at T ~ 1.4 Tc > Tstruc> which is responsible for the lattice 
instability. Our estimate of the phonon free energy, which is based on 
the method described in Ref 43, shows that the bcc lattice is favored 
over the fee lattice since the corresponding phonon free energy is 
lower by ~ 0.007 eV/ atom. This, however, is not sufficient to stabilize 
the bcc phase, since the electronic energy of the fee phase is substan- 
tially (by —0.03 eV/atom) lower than that of the bcc phase (see 
Fig. 3b). Therefore the fee lattice structure is the stable one. At even 
higher temperatures, T ~ 1.8 T^, the Tj mode becomes strongly 
anharmonic, with a pronounced double-well structure. It results in 
a high lattice entropy which lowers the phonon free energy and hence 
favors the bcc lattice structure over the fee lattice by —0.055 eV/ 
atom. This is sufficient to overcome the electronic energy loss of 
—0.04 eV/atom (see Fig. 3c) and to stabilize the bcc phase at high 
temperatures. On the basis of our results we therefore conclude that 
the high-temperature bcc (S) phase of iron is stabilized by the lattice 
entropy, which gradually increases upon heating due to the increas- 
ingly anharmonic behavior of the Tj phonon mode. 

Our results for the temperature evolution of the lattice dynamics 
offer a new perspective on the microscopic origin of polymorphism 



of iron. We find that the Tj transverse (110) phonon mode is 
responsible for the a-to-y and y-to-S allotropic transformations in 
iron. Furthermore we predict that this mode continuously softens 
upon temperature increase - a feature which is typical for simple 
nonmagnetic metals^^'"*^. For T > Tstmc ~ 1-3 Tc the bcc lattice 
becomes dynamically unstable near the N-point and collapses at 
— 1.8 Tc, at least in the harmonic approximation. By contrast, at 
temperatures T < Tstmct the phonon modes of bcc iron are relatively 
rigid, and resemble those of non-polymorphic bcc metals such as Cr, 
Mo, W. The phonon spectra of the fee (y) iron show only a rather 
weak temperature dependence^'. For the bcc phase the phonon spec- 
tra are strongly affected by temperature, i.e., the phonon frequencies 
decrease as the temperature is raised. We found that the Tj phonon 
mode softens at the fec-to-bec phase transition, which indicates the 
importance of anharmonic effects for a correct description of the 
high-temperature 5 phase of iron. Obviously, lattice anharmonicity 
and electron-electron correlations have to be taken into account 
simultaneously. This is a great challenge for future research. 

During the review process, a detailed experimental study of the 
lattice dynamical properties of iron was reported by J. Neuhaus et 
al.''^. The authors found a pronounced softening of the [110] trans- 
verse (TI) mode and conclude that the high-temperature bcc phase is 
essentially stabilized by the vibrational entropy. These new experi- 
mental findings are completely in line with our theoretical prediction. 

Methods 

We employed the LDA -t- DMFT scheme''*'"-''^ to compute lattice dynamical properties 
of iron in the temperature range covering the a-to-y and y-to-5 phase transitions. 
This state-of-the-art method combines ah initio density-functional-based techniques, 
e.g., the local density approximation (LDA) or general gradient approximation 
(GGA), with dynamical-mean field theory (DMFT) for correlated electrons systems, 
which describes the quantum dynamics of the many-electron problem exactly but 
neglects non-local effects. The LDA + DMFT approach allows one to include the 
effect of electronic correlations on the electronic band structure and phonon prop- 
erties of correlated materials. To include the electronic correlations we derived an 
effective low-energy Hamiltonian for the valence Fe sd orbitals. We employed the 
non-magnetic GGA together with plane-wave pseudo-potentials to construct a basis 
of atomic centered symmetry constrained Wannier functions for the Fe sd orbitals''*'. 
The multi-orbital Hubbard Hamiltonian obtained thereby, which includes the 
Coulomb repulsion U and Hund's exchange interaction / for the 3 d Fe orbitals, is 
solved by DMFT using quantum Monte Carlo (QMC) simulations''^. To determine 
the phase stability and lattice dynamics, we calculated the total energy and evaluated 
the corresponding phonon frequencies using the frozen-phonons method''^. The 
interaction parameters U — 1.8 eV and / — 0.9 eV are assumed to remain constant 
during the structural transformations. Further technical details about the method 
used can be found in Ref 36. 
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